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ABSTRACT 

A thermal radiative component is likely to accompany the first stages of the prompt emission of 
Gamma-ray bursts (GRB's) and X-ray flashes (XRF's). We analyze the effect of such a component 
on the observable spectrum, assuming that the observable effects are due to a dissipation process 
occurring below or near the thermal photosphere. We consider both the internal shock model and 
a 'slow heating' model as possible dissipation mechanisms. For comparable energy densities in the 
thermal and the leptonic component, the dominant emission mechanism is Compton scattering. This 
leads to a nearly flat energy spectrum {vF^ oc v^) above the thermal peak at « 10 — 100 keV and 
below 10 — 100 MeV, for a wide range of optical depths 0.03 < r^e < 100, regardless of the details 
of the dissipation mechanism or the strength of the magnetic field. At lower energies steep slopes 
are expected, while above 100 MeV the spectrum depends on the details of the dissipation process. 
For higher values of the optical depth, a Wien peak is formed at 100 keV — 1 MeV, and no higher 
energy component exists. For any value of r^ej the number of pairs produced does not exceed the 
baryon related electrons by a factor larger than a few. We conclude that dissipation near the thermal 
photosphere can naturally explain both the steep slopes observed at low energies and a flat spectrum 
above lOkeV, thus providing an alternative scenario to the optically thin synchrotron - SSC model. 
Subject headings: gamma rays: bursts — gamma rays: theory — plasmas — radiation mechanisms: 
non-thermal 



1. INTRODUCTION 

The prompt emission from gamma-ray bursts (GRB's) 
is believed to arise from the prompt dissipation of a 
substantial fraction of the bulk kinetic energy of a 
relativistic outflow, ori ginating from a central compact 
object (for reviews, see Me szaros II200I iWaxman II200I 
|Piran 2004). The dissipated energy is converted to 
energetic electrons, which produce high energy photons 
by synchrotron radiation and inverse Compton (IC) 
scattering. This model was fo und consistent with a large 
number of GRB observations ("B and et al. lllflQ^lTavanil 
Prcecc et al. 1998b). However, motivated by an 
increasing evidence for low energy spectral slopes steeper 
than the optically-thin synchrotron or sy nchrotron self 
Comp t on (SSC) in o del p re dictions llCrider et aL\ 
R997; 'Preece et al. ' '1998b^ 'Fronter a et al. I 120001 
[Ghirlanda. CcUoti & Ghisellini 2003), an additional 
thermal component was suggested which may con- 
tribu te to the observed spectrum iMeszaros fc Rees I 
20001 iMi'szaros. Ramirez-Ruiz. Rees fc Zhang I l2002t 
Rees fc Meszaros I I2005() ~ Indeed, thermal radiation 
originating from the base of the relativistic flow where 
the densities are high enough to provide (at least 
approximate) thermal equilibrium is an inevitable in- 
gredient. This radiation is advected outward as long as 
the flow (presumably, jetted) material remains opaque, 
and eventually emerges highly coUimated once the flow 
becomes optically thin, at distance Vph from the core, 
which defines the photospheric radius. The comoving 
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energy density of this thermal radiation is comparable 
to the energy density of the relativistic outflow, in the 
first stages of the flow. 

The rapid variability observed in GRB's, if inter- 
preted in the framework of the internal collisions model, 
implies that (at least some of) the internal collisions 
occur at radii sinaller than the photosphere r adius 
f Mesz aros fc Rees I l2000l iRees fc Meszaros I l2005|l . A 
similar conclusion can be drawn if alternative energy dis- 
sipation me chanisms are considered, such as magnetic re - 
connection llThomDSon 1119941 IG lannios fc Spruit II2OO5I) . 
If the energy dissipation occurs below or near the pho- 
tospheric radius, the thermal photons serve as seed pho- 
tons to IC scattering by the relativistic electrons. For 
high optical depth to scattering by the electrons (and 
the created pairs), the optically thin synchrotron-SSC 
model predictions are not valid, and a different analysis 
is needed. 

Saturated Comptonization of (synchrotron) soft 
photons was studied in the past in the con- 
text of X-rays emission from X-ray bursters 
(iSunvaev fc Titarchuk I 119801: iPqzdniakov et gP Il983t 
iZdziarski I Il98(it iZdziarski et al I Il994|) . It was first 
sugg ested as a source of GRB p rompt emiss ion by 
Liang . Kusunose. Smith fc Crider I l)1997l) . and ' Liang I 
()1997D . Comptonization by ther mal populat i on o f 
electrons was considered by Ghiselli ni fc Cellotil l)1999l ^ 
as an alternative source for GRB prompt emission. This 
model differs from the common "internal shock" model 
in that it assumes that the dissipated kinetic energy is 
continuously distributed among all the shell electrons, 
rather than being deposited into a small fraction of 
electrons that pass through the shock wave, and it 
assumes that the energy is equally distributed among 
electrons rather than following a power law distribution. 
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Th e implication of th e se mo dels were f urther analyzed 
bv iPe'er fc Waxman I 1)20041) (see also iRamirez-Ruiz I 
1200511 ^ ^ ^ 

In a previous paper ijPe'er. Meszaros k Rees 12005(1 we 
showed that clustering of the peak ene rgy at 100 keV — 
1 MeV as observed in many GRB's l|Rra,iTierd I ITflfl^ 
iPreece et al. Ill998a|l can be naturally explained as being 
due to a Compton- inverse Compton scattering balance, 
if the kinetic energy dissipation process occurs below the 
Thompson photosphere. This result was found to be in- 
dependent on the nature of the dissipation process. 

In this paper, we analyze in greater depth the effect 
of a photospheric term on the prompt emission from 
GRB's and XRF's, under different assumptions about 
the unknown values of the free parameters of the emis- 
sion models. In iJ21we show that the energy carried by the 
thermal photons is comparable to the dissipated energy 
for plausible values of the free parameters characterizing 
GRB's. This implies a large radiative efficiency, in con- 
trast to conventional internal shock models. We analyze 
in fj^lthe spectra obtained in the 'slow heating' scenario, 
and in [Qlthe spectra obtained in the 'internal shock' sce- 
nario. We show that these two different scenarios lead 
to very similar spectra for a large range of parameters 
values. In particular, we show here that a flat energy per 
decade spectrum vFi, oc v'^ arises naturally above a ther- 
mal peak at 30-100 keV. This provides a natural explana- 
tion for the ubiquitous oc prompt em ission spec- 
tra lIBand et »1 lllflM IPreece pi »1 11200(1) in the sub- 
MeV to MeV range. At lower energies, the spectrum 
depends on the details of the scenario considered; it is 
generally steeper than the typical synchrotron spectrum, 
and can become flatter through superposition of different 
regions. In [21 we present numerical results, which con- 
firm the validity of the analytical calculations in f|31 21 
and provide further details of the spectra under various 
conditions. We discuss in [JHlthe implications of the pos- 
sibility presented here of explaining the prompt emission 
by a scenario influenced by a photospheric component, 
different from the "standard" optically thin synchrotron- 
SSC model, as well as the applicability of this model to 
other compact objects, such as blazars. 

2. ENERGY DISSIPATION AND THERMAL COMPONENT 

Following the analysis of lRees fc Meszaros1l|2005ll . we 
assume that the photospheric radius rph, which is de- 
fined as the radius above which the comoving optical 
depth along the jet falls to unity, lies outside the sat- 
uration radius r^, at which the bulk Lorentz factor of 
the relativistic flow T asymptotes to the dimensionless 
entropy 77 = Lo/{AIc^), where Lq and M are the total 
energy and mass outflow rates. The saturation radius 
is rs ~ rpTy, where rp is the radius at the base of the 
flow (where the energy is released), which is compara- 
ble to the Schwarzschild radius rg of a central object of 
mass M , tq ^ avg = aUGMjc? where a > 1. Inside 
the saturation radius, the observer-frame photospheric 
luminosity is approximately the total luminosity of 
the outflow Lq, since the increasing Doppler boost can- 
cels the adiabatic decay of the comoving characteristic 
photon energy. At larger radius, r > r^, the Lorentz fac- 
tor no longer grows, and L^(r) — LQ(r/rs)^'^/'^ < Lq, the 
greater part of the energy being in kinetic form, Lk ~ Lq. 

We assume the existence of a dissipation pro- 



cess, such as magn e tic re connection l|ThomDson 1 11993: 
iGiannios fc Spr^iitl I2005D or internal shock waves 
within the expanding wind llRees fc Meszaros I 119941: 
IWoods fc Loebl 119951 ISari fc Piran I llQQTtl that dissi- 
pates a fraction of the kinetic energy. As a concrete 
example for a dissipation process, we chose the internal 
shock model of GRB's, in which variation of the flow 
Ar ~ F on a time scale At ~ tq jc results in the develop- 
ment of shocks at radius ri « 2Frs. At this radius, the 
photospheric luminosity is Lj{ri) = Lo(2F)~^/^, and the 
comoving normalized temperature of the thermal com- 
ponent is 

e = = ^ f ^0, ) ^'^ (2F)-2/3 

= 1.2 X 10-3Li/4r-5/3(a^^)-i/2^ 

where fc^, a are Bolzmann's and Stefan's constants, Lq = 
10^^L52ergs~^, F = 10^F2, and the characteristic mass 
M ~ lOmi solar masses of the central object (e.g. black 
hole) assumed, resulting in Schwarzschild radius rg = 
2GM/c^ ~ 3 X lO^^mi cm. 

The comoving proton number density in the shocked 
plasma at is Up « Lo/(47rr?cF^mpC^), therefore the 
optical depth to Thompson scattering by the baryon re- 
lated electrons at this radius is 

T^e = rronp(TT = 100L52F^5a-i. (2) 

We thus find that for the parameters values that charac- 
terizes GRB's, indeed < rph- 

We assume that a fraction of the kinetic energy is 
dissipated by the dissipation process, resulting in an in- 
ternal energy density Uint = io£ti/47rr^cF^. Fractions eg 
and £3 of this energy are carried by the electrons and 
the magnetic field, respectively, therefore the ratio of the 
thermal photon energy density Uph = aT''^ to the elec- 
tron energy density is 

A^^^ 0A4T-'^'e^X,eZ\„ (3) 

Uel 

where = lO^^e^,-!, ee = lO'^'^ee -0.5 and Uei = 
eeUint ■ Since the thermal photons provide the main pho- 
ton reservoir for electron scattering, the ratio of syn- 
chrotron and IC emitted power by the electrons is 

S = = = l.ir2' ed-ies.-o.s, (4) 

Fic Uph 

where es = 10"°-^e_B,-o.5, and ub = iBUmt- 

In order to determine the resulting spectrum, the par- 
ticle energy distribution needs to be specified. In the 
analysis below, we consider two alternative scenarios: the 
"slow heating" scenario, and the more widely used inter- 
nal shock scenario, in which power-law energy distribu- 
tion of the injected electrons is assumed. 

3. SLOW HEATING SCENARIO 

In this scenario, as suggested bv iGhisellini fc Celloti I 
1,1999,) . the dissipated energy is continuously and evenly 
distributed among the shell's electrons, during the 
comoving dynamical time of the dissipation process, 
tdyn = FAt. The electrons are assumed to have a 
Maxwellian distribution with normalized temperature 
9 el, and Comptonization is the main emission mecha- 
nism. 
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The electron temperature is determined by energy 
balance between the energy injection rate, dE+ /dt = 
Uei/rieitdyn, and energy loss rate. For optical depth 
Tje not much larger than unity, the energy loss rate is 
given by dE_/dt = (4/3)ccrT(7//3/)^itp,^(l + S). Here, 
riei is the number density of electrons in the plasma, 
nei — rip, pairs are neglected (see below), and 7//?/ is 
the electron momentum, which for a relativistic Maxwell- 
Bolzmann distribution is related to the temperature via 
^ei = 7//3j/(l + /3j)- At steady state, the electron mo- 
mentum is therefore given by 



IfPf 



1 



4T^eA(l + 5) 



1/2 



(5) 



Since A and (1 + 5) are of order unity for parameters 

— 1/2 

characterizing GRB's, 7//3/ ~ r^e • 

The resulting spectrum above the thermal peak can 
now be approximated in the following way. For 7//3/ ^ 1 
(r^e < 1), the energy of a photon after Usc scatter- 
ing is En,^ ~ (7//3/)^"='^eo, where Eq is the initial pho- 
tons' energy. In the Thompson regime, the scattering 
rate is duph/dt « riph.oneiccrT = riphflTye/tdyn, where 
nphfi is the number density of photons at the thermal 
peak. At the end of the dynamical time, the number 
density of photons that undergo Usc scattering is there- 
fore nph,nsc ~ n-phfiT^g" ■ These photons have energies in 
the range En^^--En^,^ + dEn^^, therefore 



\ ea ) ^ log (7//3/)^ j-g-j 



log I 



x(-l), 



where A not much different than 1 and S not much 
larger than 1 assumed in the last equality, resulting in 
(7//3/)^ ~ '''^e- Since Uph.n^^ is the number density of 
photons in the range En^a--En,a + dEn^^, we conclude that 
Edn/dE oc E^^, or fF^ oc i''^ above the thermal peak and 
below £max = JfiTieC^ (in the plasma frame). 

For T^e ^ 1, the electron momenta at steady state 
7//3/ ^ 1. An upper limit on the electron momenta 
7//3/ ~ 7/ is 7/ < 7 = Uei/iieimeC^ = imp/me)edeejp = 
60 erf__iee^_o.57p,0; where (7^ — l)mpC^ is the proton's 
thermal energy and jp — l7p,o (cold protons) assumed. 
This value is achieved if the optical depth is 
3 X 10"'' ^e^^-^e~^o 57~Q where A = IAq, and equa- 
tion |31 was used. For lower value of r-ye, the photons 
receive only a small fraction of the electron's energy. In 
this case only 1-2 scatterings are required to upscatter 
photons to emax = 7/TOeC^. We therefore expect a spec- 
trum with a wavy shape, composed of a low energy syn- 
chrotron peak, thermal peak, and few IC peaks (see the 
numerical results in Sj^l below). 

For high value of the optical depth, r^e larger than 
a few, the electrons momenta satisfy 7//?/ < 1. If the 
optical depth is high enough, the photons receive nearly 
all of the electrons energy, and a Wien peak is formed at 



EWP ^ WrUeC X 1 + 



6keV 



(7) 



(in the plasma frame), where equations ^ 13 were used. 
Electrons up-scatter photons up to energy (7y/3/)^meC^, 
therefore equation [T] sets a lower limit on 7//3/, 7//3/ ^ 



0.1. Equation [5I can now be used to find the value of 
r^e above which the Wien peak is formed, r^e > T^e ~ 
(3/4)(lGG/A(l 4^ S)) w 80A(7^ We thus conclude that 
for optical depth larger than 100, the spectrum ap- 
proaches a Wien spectrum, and equation provides an 
asymptotic limit of the Wien peak energy. For r^e of a 
few tens, the spectral slope below the Wien peak is softer 
than the Wien spectrum, vFy oc v'^ with a < 3. As a 
result the Wien peak is produced at energy somewhat 
higher than the asymptotic value given in equation [T] 

For values of A not much smaller than 1, equation 
implies that a large number of pairs cannot be pro- 
duced. If T~fe is larger than unity, 7//?/ < 1, the 
maximum energy of the up-scattered photons is smaller 
than TTieC^, and no pairs are produced (for synchrotron 
emitted photons to be above mgC^, a magnetic field 
B > 10^"^ G is required). For r^e smaller than 1, an 
upper limit on the number of the produced pairs can 
be obtained in the following way: For photon spectrum 
oc ly^^, assuming the upscattering photons receive all 
the electrons energy, the number density of photons en- 
ergetic than TOgC^ is approximated by n^(e > mgC^) ~ 
Uei / TTLeC^ \og['-^ f /W). The rate of pair production is 
dn±/dt ~ (3/16)ccrTn^(e > meC^)n~f{meC^ /^f S ^ S 
TTieC^), where n^{meC^ /^f ^ e ^ rrieC^) is the number 
density of photons that produce pairs with the energetic 
photons, n~f{'meC^ /^f < e < WeC^) « ^fn^^E > rrieC^) 
for Ii, oc ly^^. The number density of the produced 
pairs at the end of the dynamical time is therefore 
n± fa {Z/W)caTtdynlfuli/{meC^\og{jf/W)f, or 



riel 



16 log (i: 



6t, 



1/2 



7e 



(8) 



where characteristic value log(7//36') w 10 assumed. For 
T-ye smaller than unity this result implies that pairs can- 
not outnumber the baryon related electrons by a factor 
larger than a few at most. For 7/ < 7 = 60, a magnetic 
field larger than 10^" G is required for synchrotron 
emission at energies higher than rrieC^, thus we conclude 
that pairs are not produced by synchrotron photons in 
this case of small r^e as well. 

4. INTERNAL COLLISION SCENARIO 

We further explore the effect of a photospheric com- 
ponent on the observed spectrum in the more conven- 
tional internal shocks scenario, in which the dissipated 
energy is converted to acceleration of electrons to high 
energies. We assume a fraction epi < 1 of the electron 
number density riez is accelerated to a power law en- 
ergy distribution with power law index p in the range 
Jchar to 7max. A fraction 1 — epi of the electrons are 
assumed to have a Maxwellian distribution with normal- 
ized temperature O^i = "fchar/'2- Thermal photons are in 
the Klein-Nishina limit for Compton scattering by ener- 
getic electrons, therefore "jmax is determined by equating 
the acceleration time and the synchrotron loss time, 

7™.. = 2.3 X 10^L-V4,-V4^-l/4^^r^ay^ (9) 

where a = lao- The characteristic Lorentz factor of the 
accelerated electrons is 

eeed7p(^) 

3/2(l-ep,)+fpilog(:^) (10) 

■ 30ee,_o.5e<i,-i7p,o, 



Icha 
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where ep/ = 0.1 was taken, log(7max/7char) « 7 
and a power law index p — 2 oi the accelerated 
electrons above 'ychar is assumed. Electrons having 
Lorentz factor 7 lose their energy by Compton scat- 
tering and synchrotron emission on time scale Uoss = 
■ymeC^ /{4:/3)c(TTj^Uph{l + S), and cool down to 7/ ~ 1 
on time tioss — Uoss x 7. The ratio of the cooling time 
of electrons at ^char to the dynamical time is calculated 
using equations \213\ HI and [TUl 

tlossilchar) _ 1 _ 0-03 ^,2/3 

tdyn (4/3)-4(l + S')7c/iar-T-ye T^e ^ 

(11) 

We thus conclude that for Lorentz factor F < 2 x 10"', 
T^e > 1 directly imphes tiossilchar) < tdyn, and electrons 
lose all their energy during the dynamical time. 
If the optical depth to scattering is low, t^^ < 

0.03 Fj^^, electrons at jchar maintain their energy dur- 
ing the dynamical time. The resulting inverse Compton 
spectrum in the range OrUeC^ ■■^^har^^et? is lu on or 
vFi, = e^dn^/de oc e^. At higher energies, the spectral 
index is determined by the power law index of the in- 
jected electrons. Above ^charmeC^ Klein Nishina effect 
significantly modifies the spectrum, and numerical treat- 
ment is required (see the numerical results in below). 

2/3 

For an optical depth 0.03 Fj < r^e ^ 1, electrons 
lose their energy and accumulate at 7/ ~ 1 at the end 
of the dynamical time. The energy loss rate of electrons 
at Lorentz factor 7 < ^char is ^(7) oc 7^Wp/i(l -I- 5') oc 
7^, therefore the electron distribution below ^char is 
n-eiil) oc 7^^, and the resulting IC spectrum in the 
range OrrieC^ ..jchar'meC^ is vF^, = e^dn^jde oc with 
a = 1/2. On top of the IC component there is the 
synchrotron component from the power-law accelerated 
electrons, with spectral index a = 1 —p/2. The observed 
spectral index in this range is therefore expected to be 
a < 0.5, the exact value depends on the details of the 
scenario, i.e. the values of tpi, S and p. In ^we present 
our numerical results in this case. 

For higher optical depth, r^e > 1, electrons lose their 
energy by inverse Compton scattering and synchrotron 
emission, and accumulate at 7/ ~ 1. At low energies, 
electrons are heated by direct Compton scattering en- 
ergetic photons. We show in f0 that a photon with 
energy e > frrieC^ where / « 3 interacting with low 
energy electrons having 7 ~ 1, lose energy at an energy- 
independent rate, de/dt ~ caTnei{meC^ /2). The injec- 
tion rate of energetic photons (via IC scattering and syn- 
chrotron emission) is approximated using energy consid- 
erations, dnph{e > fmi,c'^)/dt « Uei/ fmeC^tdyn- These 
photons lose energy by downscattering to energies below 
frrieC^ on a time scale e/ (de/dt) = 2ftdyn/T^e, therefore 
at steady state the number density of energetic photons 
is given by nph(e > fnieC^) ~ 2uei/'rneC^T^e- This result 
implies that the energy gain rate of electrons via direct- 
Compton scattering is equal to the energy gain rate 
in the slo w heating mechanism, dE+ldt ~ Uei /rieitdyn 
(see also iPe'er. Meszaros fc R.eesl 120051. therefore for 
Tje not much larger than unity, equation |S1 character- 
izes the electron momenta at the end of the dynamical 
time in this scenario a well. We have previously shown 
(jPe'er. Meszaros &: Rees .2005.) that in this case of high 
optical depth, the electrons accumulate at 7//3/ — 0.1 



with a very weak dependence on the values of the un- 
known parameters, and that for optical depth higher 
than a few tens a Wien peak is formed at 1 — lOkeV. 
Equation [7| provides the asymptotic value of the Wien 
peak energy for very high optical depth in this scenario 
as well. For r^e of a few tens, the thermal peak is up- 
scattered to energy in the range Om^c? < Speak < £wp 
(see the numerical results in [jSJ. 

In contrast to the slow heating scenario, a large number 
of pairs can in principle be created, following the injec- 
tion of particles to high energies. H owever, the follow- 
ing ar gument, based on the analysis of lPe'er fc Waxman I 
f2004^ suggests that the number of pairs do not outnum- 
ber the baryon-related electrons by a factor larger than 
a few, for any value of the optical depth r^e. From en- 
ergy considerations, the pair injection rate is limited by 
dn±/dt « Uei/iTieC^tdyn- The pair annihilation rate is 

n'^caT, therefore the ratio of pairs to baryon-related 
electrons in steady state is 

^ « f . (12) 

We therefore do not expect pairs to outnumber the 
baryon-related electrons by a number larger than a few 
for any value of the optical depth. The creation of 
pairs lowers the number density of energetic photons, 
thus low ering the value of 7f /3f by a small factor (see 
IPe'er. Mfezdros fc Rees II2005D . 

5. NUMERICAL CALCULATIONS 
5.1. The numerical model 

In order to confirm the analytical calculations pre- 
sented above and to derive the spectra resulting from 
different values of the free parameters, we calculated nu- 
merically the photon and particle energy distribution un- 
der the assumptions of SI In our numerical calcu- 
lations w e used the time dependent numerical code pre- 
sented in IPe'er fc Waxman I l|2004bD . This code solves 
self-consistently the kinetic equations for and pho- 
tons describing cyclo-synchrotron emission, synchrotron 
self absorption, direct and inverse Compton scattering, 
pair production and annihilation, and the evolution of a 
high energy electro-magnetic cascade. These equations 
are being solved during the dynamical time of the dissi- 
pation process. 

In this version of the code, we assume the existence 
of a thermal photospheric component, which at the dis- 
sipation radius r.i is characterized by time independent 
luminosity = Lo{2T)~'^^^ and temperature given by 
equation ^ A fraction = O.le^,-! of the kinetic en- 
ergy is dissipated by the dissipation process. A fraction 
Ee of the dissipated energy is carried by the electrons, 
and a fraction es is carried by the magnetic field. In 
the slow heating scenario, the dissipated energy is as- 
sumed to be continuously distributed among the elec- 
trons in the two shells (and the created pairs), which 
assume a Maxwellian distribution with normalized tem- 
perature 9ei{t). The electron temperature is determined 
self-consistently at each time step by balance of energy 
injection and energy loss. In the internal shock scenario, 
electrons are assumed to be injected by the shock waves 
at a constant rate. A fraction Cp/ of the electron popula- 
tion is injected into a power-law energy distribution with 



Effect of photosphere 



5 



power-law index p between ^char and ^max (see eqs. |H1 
110(1 . The remaining fraction 1 — Cpi of the injected elec- 
trons assume a Maxwellian distribution with tempera- 
ture 9el = Iminl^ « 20ee,_o.5ed,-l • 

5.2. Slow heating scenario 

Numerical results of the observed spectra at the end of 
the dissipation process in the slow heating scenario are 
presented in Figure ^ for different values of the optical 
depth T-ye- In producing the plots in this graph, we chose 
representative values of the free parameters, a = 1, T = 
10^, e<j = 10-1, ee = es = lO^^ ^ The differences in 
the plots are due to different intrinsic luminosity, Lq = 
IQ-is _ iq53 ergs^i which result in a wide range of optical 
depths, 10~^ < < 10'^ (see eq. EJ- A similar graph is 
obtained by assuming different value of T and constant 
luminosity, where variation in the value of F by a factor 
of 3, leads to variation in r^e by a factor 10^'^. 

In the case of low optical depth r < 1, one expects the 
thermal peak at the observed energy 30FmeC^/(l -f z) « 
10 — 100 keV, caused by the adverted thermal photons. In 
the scenario of r = 0.01, the characteristic Lorentz fac- 
tor of the electrons is jf w 10, and synchrotron radiation 
produces the peak observed at low energy, ^ 300 eV. For 
higher optical depth the characteristic electron Lorentz 
factor is lower, 7/ > 1 (see eq. 13)), and the synchrotron 
emission occurs at lower energies, below the synchrotron 
self absorption frequency (see the analysis of synchrotron 
self-absorption frequency as a function of the free param- 
eters values in lPe'er fc Waxman 11200'^ . 

The high value of the characteristic Lorentz factor of 
the electrons in the scenario with r = 0.01 implies that 
two scatterings are required to upscatter thermal photons 
to ^ 100 MeV, therefore Compton scattering of the syn- 
chrotron and thermal photons produces the wavy shape 
of the spectrum above 10 kcV . 

At moderate values of the optical depth, r ~ 0.1 — 10, 
multiple Compton scattering produces the flat spectrum 
observed above the thermal peak, in the range 10 keV — 
10 MeV, in agreement with the predictions of ^ At even 
higher values of the optical depth, r-^e ^ 100, a Wien 
peak is formed at sub MeV energy. As explained in 
iPe'er fc Waxmani 1(20041) . the high optical depth imphes 
that in this case photons are not observed prior to a pe- 
riod of adiabatic expansion, during which 50-70% of the 
photons' energy is converted to kinetic energy (bulk mo- 
tion) of the plasma. The observ ed Wien peak is therefore 
expec ted at « 300 keV (see also [Pe'er. Meszaros fc Reesi 

5.3. Comparison of scenarios: low optical depth 

Figure El shows comparison of the spectra obtained 
for the two dissipation scenarios considered (slow heat- 
ing and internal shocks) and two different values of the 
magnetic field, in the scenario of very low optical depth, 
Tje = 4 X 10"'^. For this value of the optical depth, only 
a small fraction of the electrons energy is radiated. If the 
magnetic field is weak (dash-dotted lines), the adverted 
thermal peak at 10 keV is prominent. In the slow heating 
scenario (thin dash-dotted line) an IC peak at ~ 15 MeV, 
created by electrons having 7//3/ — 30 is expected. The 
thin solid line, representing the slow heating scenario 
with equipartition magnetic field, shows four peaks: The 
synchrotron peak at 300 eV, the thermal peak at 10 keV 



and two peaks produced by IC scattering the low energy 
peaks, at 300 keV and 10 MeV. The combined effect of 
these peaks is a wavy shaped, nearly flat spectrum in the 
range 100 eV - 100 MeV. 

An internal shock scenario with weak magnetic field is 
shown by the thick dash-dotted line in Figured Elec- 
trons injected with a power law index p = 2 between 
jchar and jmax producc the high energy component 
by IC scattering the thermal photons. Electrons with 
Lorentz factor 7 < 7 = (46') = 3 x 10^ cool by IC 
emission faster than the dynamical time, therefore the re- 
sulting IC spectrum is nearly flat, i/F^, = e^dn^/de oc e" 
with a = below Tj'^emeC^ /{I + z) 30GeV. For en- 
ergetic electrons with 7 > 7, the thermal photons are in 
the Klein-Nishina regime, and as a result these electrons 
cool slower than the dynamical time, and the spectral 
slope of the IC spectrum at higher energies is a = 1/2. 
If magnetic field is added (thick, solid curve), electrons 
at all energies cool fast by synchrotron emission, and 
the resulting synchrotron spectrum is flat in the range 
lOOcV- lOGeV. 

5.4. Comparison of scenarios: intermediate and high 
optical depth 

For low to moderate values of r^e, equation [S] provides 
a good estimate of the electron momenta. The ratio be- 
tween the energy of the synchrotron peak Ssyn and the 
energy of the thermal peak eth = iOrUeC^ is given by 

Ssyn _ 3 hqB {-ffPfY 
eth 2 meC Eth 

^ 21/3 Hq rV3(eae,)i/2e,^ 

V 8 / k% meC A{1 + S)Tje 

_4^2<l(r^ pl/3 1/2 1/2 .-1 ofc. /.o^ 

where = lO^'e"^ 4 eV is the observed energy of the 
thermal peak. We therefore conclude that for t^^ ^0.1 
the peak energy of the synchrotron component is much 
lower than the peak energy of the thermal component, 
and is below the self absorption frequency for Tje ^ 1- 
This conclusion is confirmed by the numerical results of 
the slow heating scenario (thin line) with strong (solid) 
magnetic field, presented in Figures |21 and 01 for r-y^ — 
0.1, 1 respectively. A low energy synchrotron component 
at 200 eV is prominent for r^e = 0.1 (Figure^, and is 
much weaker for x^e = 1 (Figure Ej). 

Figures I3I5I provide a comparison of the spectra ob- 
tained for the two dissipation mechanisms, equiparti- 
tion vs. low magnetic field, and two different values 
of epi. The internal shock scenario is shown with the 
thick lines in Figure Ol The solid line denote the sce- 
nario of equipartition magnetic field and epi = 0, i.e., 
the electrons are injected at jchar and than cool to 7/. 
Synchrotron emission dominates the spectrum below the 
thermal peak, at 0.1 — lOkcV, and IC scattering domi- 
nates at higher energies. The electron distribution below 
Ichar is riei{'^) cc 7~^, therefore the emitted radiation 
has spectral slope a = 1/2. IC scattering by electrons 
accumulated at 7/ results in somewhat flatter (a < 1/2) 
spectral slope. The thick dash-dotted curve shows a sce- 
nario with weak magnetic field (e_B = 10~^) and epi = 1. 
In this scenario, energetic electrons upscatter photons to 
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high energies; however, pair production suppresses emis- 
sion above GeV. The number of pairs created in this 
scenario is found numericaUy to be rt±/nei ~ 4, in agree- 
ment with the analytical predictions of 21 

The numerical results obtained for r-ye — 1, 10 are pre- 
sented in Figures 01 and These results further show 
the similarity between the spectra obtained by the two 
different dissipation mechanisms, as well as the weak 
dependence of the spectral shape on the value of the 
magnetic field or on e^/. Multiple scattering creates the 
flat energy per decade spectrum above the thermal peak 
at WYmeC^I[\ -f z) ~ 30 - 100 keV, for a typical T = 
100 — 300. This energy is very similar, albeit somewhat 
lower than the observed spectral break energy, ~ 300 keV 
(|Band et al. liTooliPreece et al. Ill998all2000|) . The flat 
spectrum extends up to T{jfPf)mec'^ /{I -I- z) ~ 10 MeV. 

In the internal shock scenario, the energetic photons 
produce pairs, resulting in a cutoff above ^ 100 MeV. 
The number density of pairs in the scenario r^e = 1 
is found numerically to be n±/nei — 4 (ep/ = 0), 
n±/nei ~ 5 (cp/ = 1). Similar, but somewhat lower 
values [n±/n(,i — 1.5,3.5, respectively) are obtained 
for T^e = 10. We thus find that, indeed, pairs do 
not outnumber the baryon related electrons by a large 
factor, due to pair annihilation. The pair annihila- 
tion phenomenon causes t he small peak observed at 
rmecV(l+z) = 25 MeV fsee lPe'er fc Wa.xma,n 120011 for 
further discussion). In summary, for intermediate values 
of the optical depth r^g w 0.1 — 10, multiple IC scat- 
tering produces an approximately flat energy per decade 
spectrum in the range « keV - sub GeV for the two dis- 
sipation mechanism considered regardless of the values of 
the magnetic fleld or of ep/. 

Figure El shows the asymptotic case of high opti- 
cal depth, T^e = 100, 1000. For optical depth larger 
than ~ 100, multiple Compton scattering leads to the 
formation of a Wien peak at sub- MeV energies, re- 
gardless of the details of the acceleration mechanism, 
or the values of any of the free parameters (see also 
iPe'er. Meszaros fc R.ees l)2005^ . As shown in 21 for high 
values of the optical depth the number density of pairs 
is strongly suppressed by pair annihilation. This calcu- 
lation is confirmed by the numerical results, n±/nei — 1. 

6. SUMMARY AND DISCUSSION 

In this work, we have considered the effect of a pho- 
tospheric component on the observed spectrum after ki- 
netic energy dissipation that occurs below or near the 
thermal photosphere of the outflows in GRB's or XRF's, 
with particular emphasis on the resulting spectral shape. 
Two dissipation mechanisms were considered: the "slow 
heating" mechanism, and the internal shock scenario. 
We showed that the resulting spectra are largely inde- 
pendent on the details of the dissipation mechanism or 
of the values of any of the free parameters. A strong 
dependence was found only on the value of the optical 
depth, T^e- For r-ye smaller than a few tens, an approxi- 
mately flat energy per decade spectrum is obtained above 
a break energy at few tens - hundred keV and up to sub 
GeV, while steep slopes are obtained below the break en- 
ergy (see Figures 01 ISJ . For higher values of the optical 
depth, a Wien peak is formed at sub MeV energies, with 
steep slope at lower energies, and sharp cutoff above this 
energy. The production of pairs does not change this re- 
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Fig. 1. — Time averaged spectra obtained for different values 
of optical depth, r^e in the slow heating scenario. Results are 
for a = 1, r = 100,ed = 0.1, Ce = IQ-'^-S and GRB luminos- 
ity L = 10-"', 10-3, 10-2, 10-\ 10°, 10\ 10^, with corresponding 
optical depth r^e = 10-^ - 10^ (see eq. 0. <lb = 10~° '' for 
L = 10-* - 10°, and eg = 10-^ for L = 10^,102. A redshift z = 1 
in a flat universe with Qrn = 0.3, Ha = 0.7, Hq = 70 is assumed. 
The spectrum is not corrected for the energy loss due to adiabatic 
expansion following the dissipation, which lowers the energy for 
the cases of L = 10^, 10^ (r-ye = 10^, 10^) by a factor of 2-3. 




Obs. energy [eV] 

Fig. 2. — Time averaged spectra obtained for low value of the 
optical depth, t^^ = 4 X 10-^. a = 1, L = 10-^, F = 300, 
e^i = 0.1, ee = 10— "'^ assumed. Results for slow dissipation are 
shown in thin lines and dissipation by shock waves with power law 
index p = 2 and Epi = 1 are presented in thick lines. Solid lines 
arc for high magnetic field, eg = 10-°'^ and dash-dotted lines are 
for es = lO-''. 2 = 1 with the same cosmological parameters as in 
figure m 



suit, since pair annihilation limits the ratio of the number 
density of pairs to baryon-related electrons to a factor not 
larger than a few. 

Our theoretical and numerical results emphasize the 
important role of inverse Compton scattering in the 
formation of the spectra in scenarios involving an im- 
portant thermal component. As a consequence, the 
observed spectral slopes in the two model consid- 
ered are significantly different than the synchrotron 
model prediction. The (nearly) thermal component 
observed for low to intermediate values of the opti- 
cal depth, and the Wien peak formed for high value 
of the optical depth can account for the increasing 
evidence for steep spectral slopes observed below ~ 
100 keV in the early stages of GRB's l|Preece et al. 120021: 
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Fig. 3. — Time averaged spectra obtained for low value of the 
optical depth, r^^ = 0.1. a = 1, L = 10"^, T = 100, = 0.1, 
ee = 10""'^ assumed. Results for slow dissipation are shown in thin 
lines and dissipation by shock waves with power law index p = 2 
are presented in thick lines. Solid lines are for high magnetic field, 
(.B = lO"'''^ and dash-dotted lines are for eg = 10 . £p; = 
(thick, solid), and e^; = 1 (thick, dash-dotted). 2 = 1 with the 
same cosmological parameters as in figure ITl 
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Fig. 5. — Time averaged spectra obtained for intermediate value 
of the optical depth, r-ye = 10. Results are shown for L = 10~^ 
with all other parameters same as in Fig. |3] Results for slow 
dissipation are shown in thin lines and dissipation by shock waves 
with power law index p = 2 are presented in thick lines. Thick, 
solid line: high magnetic field, es = 10""'^ and epi = 0. Thick, 
dash-dotted line: eg = 10~^, tpi = 1. Thin, dash-dotted: slow 
heating, eg = 10~®. 
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Fig. 4. — Time averaged spectra obtained for low value of the 
optical depth, r^e = 1. Results are shown for L = 10~^ with all 
other parameters same as in Fig. |21 Results for slow dissipation 
are shown in thin lines and dissipation by shock waves with power 
law index p = 2 are presented in thick lines. Solid lines are for 
high magnetic field, eg = 10"*^'^ and dash-dotted lines are for 
eg = 10~®. epi = 1 (thick, solid), and epi = (thick, dash-dotted). 
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Fig. 6. — Time averaged spectra obtained for high value of the 
optical depth, r-^e = 100, 1000. Results are shown for L = 1 (solid), 
L = 10 (dash-dash) with all other parameters same as in Fig. 121 
Results for slow dissipation are shown in thin lines and dissipation 
by shock waves with power law index p = 2 are presented in thick 
lines. Thick, solid line: high magnetic field, eg = lO"'''^ and epi = 
1. Thick, dash-dash line: eg = 10~®, epi = 0.1. Thin, solid: slow 
heating, eg = 10""'^. Thin, dash-dash: slow heating, eg = lO"''. 



IGhirlanda Celloti fc CAise lhrn I llffll IRvde 1120041) . As 
shown in iPe'er. M( 3szaros fc R,ees I II200I " and as con- 
firmed here (figures H El ; for a scattering optical depth 
larger than ^ 100, dissipation at a sub-photospheric ra- 
dius can naturally lead to a clustering of the peak en- 
ergy at ^ 100 keV — 1 MeV, as observed in many bursts 
ijBrainerd 111991 iPreece et al. Ill998all2000|) . 

For low values of the optical depth, r^e ~ 0.1 — 1, 
the advected thermal component produces an observed 
peak at a few tens of keV. This peak energy is con- 
sistent with t he peak ener gy observed in X-ray flashes 
at ~ 25keV l|Heise et al. m OOl') and X-ray-rich GRB's 
(XRR's) , which show p eak energy clustering at ~ 50 keV 
((Sakamoto et aZ. ]|2005|) . Thermal peak energies in this 
range of few tens of keV are obtained for a wide range 
of values of the free model parameters, lO^'' < L52 ^ 1, 
1 r2 < 3 and 1 ;$ a ;$ 10^ which result in an op- 



tical depth ~ 1 (see eas lll2ll . It is difficult to obtain 
in our model a thermal peak energy at energies below 
a few keV combined with an optical depth larger than 
^ 10~^. Therefore, for a peak below a few keV, elec- 
trons maintain their energy (see 21 . We conclude that 
if XRF's peak energies are due to advected thermal com- 
ponent, no XRF's with peak energy below a few keV are 
expected. 

The inclusion of an advected thermal component in 
the prompt emission calculations can help resolving the 
question of the accelerated electron distribution, which 
is of high theoretical importance. As pointed out by 
faring & Brabv ( 2004), fitting the data of many GRB's 
by synchrotron and Compton emission only, require the 
acceleration of the majority of electrons to a power law 
distribution, leaving far too few thermal particles than 
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to be expected as a seed distribution for the power law 
population. Here, we overcome this problem by assuming 
that the th ermal component seen in many bursts (e.g., 
iRvde 1120041) . and presumably exist in all of them at early 
times, is advected from the core. As a result, a vari- 
ety of spectra can be obtained under different assump- 
tions on the values of the free parameters, without the 
requirement that most of the electrons are accelerated to 
a power law distribution (in the presented figures of the 
power law acceleration model, we allow the parameter 
tpi to vary between and 1). 

The spectra presented in this paper ([JSJ describe the 
emission resulting from a single dissipation phase (e.g., 
single collision between two shells) for a particular choice 
of model parameters. Observed spectra are expected to 
be combination of spectra produced by many such dis- 
sipation processes, which are characterized by different 
parameters, e.g., different locations of collisions within 
the expanding wind. A detailed comparison with obser- 
vations therefore requires a detailed model describing 
the distribution of single dissipation par ameters within 
the fireball wind (see, e.g., [P aignc fc Mochkoyitch I 
T99^ IPanaitescu. SnadTfc Meszar^ 

Guetta. Soada fc Waxmanl |200H) . The construc- 



tion and investigation of such detailed models would 
involve additional assumptions and parameters, which 
is beyond the scope of this manuscript. 

Nonetheless, the similarities found between the re- 
sulting spectra, produced under reasonable assumptions 
about the unknown values of the free parameters, sug- 
gests that the key results would remain valid for a spec- 
trum that is obtained from a series of dissipation pro- 
cesses. These results include an observed break energy 
at few tens of keV to sub MeV, accompanied by slopes 
steeper than the optically thin synchrotron spectrum at 
lower energies, which are obtained for a single dissipation 
event for any value of the optical depth, and flat energy 
per decade spectra obtained up to sub-GeV energy for 
optical depths of 10^^ - few tens. For a superposition of 
dissipation events at various radii, flatter slopes can be 
obtained. 

A superposition of dissipation events is also required 



in our model in order to explain both peak clustering at 
~ 300 keV and flat energy per decade spectra at higher 
energies, as observed in many GRB's. A high optical 
depth r^e ^ 100 is required in order to obtain a Wien 
peak at sub-MeV energies, while lower values of r^e are 
required in order to obtain a flat ener gy per decade spec- 
tra at higher energies. As detailed in ijPe'er fc WaxmanI 
l2004HPe'er. Meszaros fc Rees 120051) . for the case of high 
optical depths the photons are expected to lose ~ 30% of 
their energy to the bulk motion of the plasma during an 
adiabatic expansion phase before escaping, thus the ob- 
served Wien peak energy is expected at ~ 300 keV. We 
conclude that given such a superposition, a clustering of 
the peak energy, steep slopes below the peak and flat 
energy per decade spectra above the break are a natu- 
ral consequence of a model containing a significant ther- 
mal emission component. These results show good agree- 
ment with a large number of observations of GRB prompt 
emission sp ectra l|Band et al. 11199.4 iPreece et a,l. 1120001 
IRvde I l2007). 

The fluxes predicted by the model are within the detec- 
tion capability of the Swift^ satellite in the keV range, 
and the GLAST^ satellite in the sub- GeV range. Obser- 
vations of a cutoff at high energies will therefore provide 
information about the optical depth during the dissipa- 
tion phase, hence constraining the value of F, one of the 
least constrained free parameters of the model. 

The results may be applicable to a variety of com- 
pact objects, such as GRB's, XRF's and blazars. At 
least 60 flat spectrum radio quasars having flat (a w 
0) spectral index in the EGRET range were reported 
(e.g.. iMukheriee et al. Hl997j) . Similar spectra were ob- 
tained for several quasars using the BeppoSAX satellite 
ijTavecchio 1120001) . while clusteri ng of blazars peak ener- 
gies at 1-5 MeV was reported bvlMcNaron-Br own et aL\ 
( 1995). We thus conclude that similar mechanisms to the 
ones presented here may explain some of the observable 
effects in radio quasars as well. 



AP wishes to thank R.A.M.J. Wijers for useful discus- 
sions. 



http://www.swift.psu.edu 
^ http://www-glast.stanford.edu 



APPENDIX 

ENERGY LOSS RATE OF HIGH ENERGY PHOTONS DUE TO COMPTON SCATTERING WITH 

MILDLY RELATIVISTIC ELECTRONS 

The rate of scattering by a single electron having Lorentz factor 7 passing through space filled with a unit den- 
sity, isotropica lly distributed, mono-ener getic photons with energy aivrieC^ was first derived bv ..Tones (1968) . and 
summarized in iPe'er fc WaxmanI l)2004bl) , 



dtda 



[F{C+) - F(C- 



(Al) 



2U/2 



Here, a is the energy of the outgoing photon in units of rrieC , tq is the classical electron radius, /3 = (1 — I/7 ) 
(± are the upper and lower integration hmits and F( C) is given by the sum of 1 2 functions, F{() = J2]=i /i(C)j where 
the functions /i(C) are summarized in, e.g., eq. 20 of lPe'er fc WaxmaiTI lj2004b|) . 

The limits C±(a;7iQ;i) depend on the energy of the outgoing photon, a for given ai and 7 f see IPe'er fc WaxmanI 
|2004b). For ai ^ 1 and 7 — 1, this dependence is degenerated and (± = 1 ± (3. In this approximation, in calculating 
the difference of the 12 function Afi = /i(C+) — /i(C-) only 4 out of the 12 terms are to the order aj~^ or a~^, while 
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Fig. A7. — Energy loss rate of photon with initial energy aimeC^ due to Compton scattering with unit density, isotropically distributed, 
mono-energetic electrons with velocity /3 = 0.1. 



the other 8 are to lower orders in a, ai. Thus, to leading order in a, ai, 

A/3^/3(C+)-/3(C~) ^ 



2/?(l+r) 
, a ' 



(l-/3^)ai' 



(A2) 



A/i 



2/3(l+/3") 
771 



ai 



1 



1 



The (normalized) energy loss rate of photon being scattered by unit density, isotropically distributed, mono-energetic 
electrons at Lorentz factor 7 is therefore given by 



da{-f,ai) 
dt 



Y 

dtda' 



da' 



cctt 



1-/52 



3^ 
16 



CCTt 



5^ 
16 



(A3) 



In order to check the validity of this approximation for values of ai not much larger than 1, a numerical integration 
of eauation lAll was carried out, using the exact formula. The results of the energy loss rate are presented in Figure 1X71 
From the figure we find that the approximation de/dt « caTTdi^c? where e = arrieC^ for the energy loss rate of 
photon due to Compton scattering with unit density, isotropically distributed, mono-energetic electrons with 7 « 1 is 
valid for e > frrieC^, with / « 3. 
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